Coulomb analogy for nonhermitian degeneracies near quantum phase transitions 
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Degeneracies near the real axis in a complex-extended parameter space of a hermitian Hamiltonian 
are studied. We present a method to measure distributions of such degeneracies on the Riemann 
sheet of a selected level and apply it in classification of quantum phase transitions. The degeneracies 
are shown to behave similarly as complex zeros of a partition function. 
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Quantum phase transitions (QPTs) appear in infinite 
lattice systems (like spin arrays with bound range of in- 
teractions) [l| as well as in finite many-body systems of 
interactin g bo sons or fermions, see e.g. Refs.0, H, 0, H, 
H S H, i,M [O, [H, [il- QPTs are driven by interac- 
tion parameters of the Hamiltonian and may affect the 
ground state 3, [ p , B, H . H, H, US] as well as individual 
excited states [ll|, llJ, Although the QPT nonana- 
lytic features strictly occur only in the system's infinite 
size limit, significant precursors can be observed already 
at a moderate size. 

Typical examples of finite many-body systems with 
QPTs are provided by two-level bosonic models, where 
one of the bosons, s, is scalar, while the other, 6^^^ 
transforms under rotations as the L'^'rank tensor p^ . 
The single-particle Hilbert space has a dimension 2L+2, 
hence the term "finite" . Such models with i = 0, 1, 2 are 
used to describe various types of collective excitations in 
atomic nuclei and molecules [l5| . In absence of in- 
teractions between s and b bosons, the ground state is a 
pure one-component condensate, usually the one with s 
bosons (the system is in its "spherical phase"). As the 
interaction strength ^ increases, a certain critical point 

is reached where the ground state becomes a mixture 
of both types of bosons (the "deformed phase" arises) . 

An important distinction between the infinite and fi- 
nite models follows from the fact that in the latter ones 
a value oc A^"^, where TV is the number of particles, can 
be identified with the Planck constant. Hence iV— i-oo is 
the classical limit [2]. Although this implies that QPTs 
in finite many-body models are rooted in semiclassical 
properties, there exists a plethora of genuinely quantum 
signatures. Some of them will be discussed in this Let- 
ter, where a link is elaborated between the occurrence of 
various types of QPTs and the distribution of degenera- 
cies of Hamiltonian eigenvalues in the complex-extended 
parameter space. 

We will consider a class of Hamiltonians depending 
linearly on a single dimensionless control parameter A, 



H{X) = Ho + XV, 



(1) 



where [iJo,^^] 7^ 0. The Hilbert space dimension n is 
finite and energy eigenvalues Ek (A) are counted from k=0 



(the ground state) to A: = n— 1. As generally known, the 
maximal rate of change of a given eigenvector |i/'fe(A)) 
appears when the level k undergoes an avoided crossing. 
A QPT takes place if the relative spacing between levels 
(in units of the mean spacing) becomes infinitely small 
as iV — > cx). If the limiting process results in a crossing 
of two levels, the QPT is of the first order (with a jump 
of -^Ek and the corresponding swap of wave functions). 
On the other hand, if the number of interacting levels is 
locally large but with no real crossing, the resulting QPT 
is continuous (with a softer type of nonanalyticity) . 

The occurrence of avoided crossings in the spectrum 
signals the presence of true degeneracies somewhere 
nearby, in the plane of complex-extended parameter A = 
X+ifj, 0, [lB|- The generalized Hamiltonian H{A) is non- 
hermitian and requires to handle the left and right eigen- 
vectors. The eigenvalues live on n Riemann sheets which 
can be labeled by the ordinal number k of the respec- 
tive level on the real axis. There can be two simplest 
types of degeneracies: (i) a diabolic point, when the two 
sheets just touch each other [l3], and (n) an exceptional, 
or branch point [l^, if the two sheets are entangled by 
the square-root type of singularity [l^, [13, [2l|. In our 
case, the common type of degeneracy is (ii). This fol- 
lows from a perturbative expansion near (but not at) the 
degeneracy which allows to approximate the evolution of 
the two close eigenvalues by a 2x2 matrix with generi- 
cally nonzero offdiagonal elements. In this situation, the 
degeneracy is a branch point. A diabolic point would 
require an additional constraint that offdiagonal matrix 
elements vanish at the singularity 21]. Since in the fol- 
lowing we do not check individual nature of each crossing, 
we use a general term "nonhermitian degeneracy" [2^ . 

If a nonhermitian degeneracy comes close to the real 
A axis, an avoided crossing is encountered in level dy- 
namics of H{X). We will therefore focus on the distri- 
bution of degeneracies in the A plane. The degeneracies 
are simultaneous roots of the characteristic polynomial 
P{E,k) = dci[E - H{K)\ and its derivative -§gP{E,K). 
The elimination of variable E leads to the discriminant 
D{K) (proportional to the resultant), a polynomial of or- 
der rt(n— 1), whose roots indicate positions Aq. of the 
degeneracies. For real Hq and V , the roots come as 
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TVd = hn{n — 1) complex-conjugate pairs, so 



2?(A)(x n(A-A„)(A™A:) 



(2) 



The discriminant can be also expressed as [11 



D{A) l[[E,{A) E,{Ar = {-r-l[D,{A) , (3) 

k<l k 

Dk{A)^ l[[Ei{A)-Ek{A)], (4) 

where we introduced a factorization by "partial discrimi- 
nants" Dk ■ The (unknown) functions Ei (A) denote com- 
plex eigenvalues on individual Riemann sheets, which are 
pairwise connected by degeneracy points Aki—Aik defined 
by Ek{Aki)=Ei{Aki). The set {Am} with fc</ is equiva- 
lent to {Aq} with a—1, . . . , A^D- 

If analyzing the /c*'^state avoided crossings, only the de- 
generacies on the fc'^ Riemann sheet are relevant. These 
are zeros of D^.. However, with an increasing dimension 
n the assignment of degeneracies to individual sheets be- 
comes practically impossible. Fortunately, there is a way 
how to indirectly measure the distribution of fc'^sheet de- 
generacies close to the real A axis. Squaring Eqs. (HJ and 
^ we find that is a polynomial of order 2(n— 1) with 
the roots A^; and At,. On the real axis, A = X + iO, one 



^kl 

therefore has 



DkiX)' 



^kl 

i(^k)^ — - — 



(5) 



Ck > 0, where Rki{X) denotes the distance of the point 
Aki = Xki + ifJ-ki from place A. Hence Dk{XY is sensi- 
tive to the proximity of the fc^'^sheet degeneracies to the 
real axis and, at the same time, it is expressible through 
Eq. (IH), using solely real eigenvalues of Hamiltonian ([T]). 
By defining Uk{X) = — ^ InDfc(A)^, where O = n— 1 is a 
convenient scaling constant, we replace the product over 
degeneracies by a sum: 

ii=ik) 

In the first line we identify a 2-dimensional Coulomb po- 
tential (scaled and shifted) caused by n— 1 point charges 
at positions {Xku^^ki) in the complex plane. Therefore, 
the introduction of Uk makes it possible to benefit from 
the basic intuition of planar electrostatics. (Note that 
energy levels in the second line can also be thought as 
point charges.) 

If the fc*'^eigenstate exhibits a QPT at a certain critical 
point AJ: on the real axis, the degeneracies are expected 
to approach infinitely close to this point in the N ^ oo 
limit 0, . Since the QPT depends most sensitively on 



the closest degeneracies, it is convenient to define the 
following additional quantities: 

Fk{X)^~^UkiX), Ck{X) = ^Fk{X) , (8) 
gfe(A)=^hm|^ Ck{X')dX' = \un^[Fk]ltl- (9) 



Fk is nothing but a repulsive force felt by a "trial charge" 
on place A due to the charges in the complex plane. 
As the charge distribution is symmetric against complex 
conjugation, the force is (anti)parallel with the real axis. 
If the charges approach very close to A^,, the force will 
change its sign from — to -I- at this point and Ck, which 
represents a change rate of the force, will have a max- 
imum there. Sometimes the force may even jump at 
A^, i.e. to fiip its sign without actually crossing zero. 
Such a place corresponds to asymptotic "condensation" 
of charges on the real axis (a locally 1-dimensional charge 
layer across ImA = is formed) and can be detected by a 
nonzero value of Qk in Eq. ([9]). In contrast, Qk is obvi- 
ously zero as far as the force is continuous. 

It turns out that the quantities defined in Eqs. ([S])- 
^ are useful for the classification of the fc^'^level QPTs. 
From Eq. (O we get 
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4>kl 



A.[EiiX) ~ EkjX)] 
\EiiX) ~ Ek{X)\ 



(10) 



where 4>ki = for Ei > Ek and (pki = 1 otherwise. This is 
a continuous function of A unless some of the derivatives 
^ [El — Ek] are discontinuous. Since the discontinuity 
(if any) appears in the N—i-oo limit, when the spectrum 
becomes infinite, a "macroscopic" fraction of nonanalytic 
energy differences needs to sum up to produce a finite 
jump of Fk at a certain point. This is naturally satisfied if 
j);Ek itself has a jump, i.e. if the fc^'^state exhibits a first- 
order QPT. Therefore, Qk 7^ for a generic first-order 
transition while Qk — for a continuous transition or for 
a crossover. The latter two cases can be distinguished 
by the behavior of Ck, which shows an asymptotically 
singular peak at A^ if there is a QPT of any type. 

To illustrate these matters, let us give two examples. 
In the first one we assume that the A:*'^sheet degeneracies 
are all located on a line perpendicular to the the real A 
axis at A = A^ and that they become infinitely dense as 
N—KX). Introducing (5 = A — A^ we obtain 



Cki6) = / pfe(M) 



-52 



(11) 



where pk {p) is an asymptotic linear density of degenera- 
cies normalized to unity: /9fc(/i) = limjv^oo ji'Ylii^it^~l^ki)- 
Note that the normalization is ensured by the factor l/fi 
which takes the reciprocal value of the number of de- 
generacies. In any case, Ck has a symmetric maximum 
at (5=0 where all its odd derivatives (if they exist) van- 
ish. It becomes nonanalytic if the support of /Ofc(^) has 
its infimum at /i=0. In particular, let us assume that 
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FIG. 1: Ck = -iirt^fe, see Eq. (0, obtained from the J^v^O 
spectrum of Hamiltonian (|12p for some selected values of fc = 
xn and A'^. The peaks indicate the passage of a given level 
through an avoided crossing at E=Q [23{ . 
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close to this point the density can be approximated by 
Pkip) ^ fJ-^^ with p>0. Then we can show that: Qk{0) 
is finite for p — and zero otherwise, Cfc(O) diverges 

for p e(0,l] and is finite for p >1, ^Cfe(O) diverges 

for p€{l,3] and is finite for p>3, ^Cfc(O) diverges for 
pG(3,7] and is finite for p>7 etc. Therefore, the line ge- 
ometry of degeneracies induces a first-order QPT if p=0, 
i.e. lim^^o Pfe(A') 7^ 0. In all the other cases the QPT is 
of a continuous type. An increasing power p shifts the 
6=0 nonanalyticity to higher and higher derivatives of i*!- 
(and Ek). This sheds light on the asymptotic nature of 
level crossings in the first-order and continuous QPTs. 

The second example shows the behavior of nonhermi- 
tian degeneracies in the interacting boson model of nu- 



clear physics jl4| . It is one of the above-mentioned sb 



bosonic models in which b^^^ is identified with the L — 2 
boson called d. The ground-state QPTs between spheri- 
cal and deformed phases can be of both first and second 
order ^. Here, we will focus on the [0(6)-U(5)]dO(5) 
transition [2^, which exhibits a second-order QPT for 
the ground-state and a chain of continuous QPTs for ex- 
cited states [H, [l3| . 

The Hamiltonian per boson reads as 



+2 



(A) = A inrf - (1 - A) ^ }^ H"'QmQ-m , (12) 



where Ud stands for the d-boson number operator and 
Q„ = dl^s + (-)'"std_^. Eq. IH]) has the form (P) 
with A G [0, 1], the value ^ = 1 — A rising the strength of 
interactions between s and d bosons. The ground-state 
QPT is located at Aq = |, where the classical poten- 



tial corresponding to Hamiltonian changes between 
sombrero-like (A < Aq) and quartic oscillator (A > Ag) 
forms . Excited-state QPTs take place for A < Aq at ab- 
solute energy £'=OjWhere the sombrero potential has the 
central maximum [12|. As individual levels cross the top 



FIG. 2: Large- A'' evolution of the a;=0 and 2;=0.1 peaks of 
Ca:(A) for Hamiltonian (|12|) . Log- log plots from top to bottom 
show the peak height (relative units), inverse width at half 
maximum, and a product height x width. Lines indicate a 
power-law asymptotics for 7V> 10*. 



of this maximum, wave functions become singular and 
level energies evolve with locally infinite curvatures p^ . 
This is related to the crossing of the classical phase-space 
separatrix, which exists only in absence of the centrifugal 
barrier, i.e. for states with vanishing 0(3) and 0(5) rel- 
ative "angular momenta" J/Jmax and u/wmax [131 ■ For 
finite N, the J=v=0 level dynamics shows a character- 
istic pattern of avoided crossings propagating through 
the spectrum at £'=0 23]. The passage of a given level 
k through the crossing yields a peak of Cfc(A), as illus- 
trated for some selected levels in Fig. [T] Here we replace 
integer k by an excitation ratio x — k/n £ [0, 1]. The 
peak centroid moves leftwards linearly with increasing x 
and sharpens with an increasing boson number N. 
From semiclassical arguments it is known that the 

QPTs of Hamiltonian are continuous, with discon- 
,2 

tinuous or infinite -j^Ek for a:=0 or x >0, respectively 
[13, [3 ■ In finite- calculations, however, one can never 
verify whether is zero at A^ or not. To address 
this problem, we determined the dependence Cx{X) for 
x=0 (the ground state) and a;=0.1 (excited state in 10% 
of the spectrum) for boson numbers up to 10^. From 
these calculations we extracted the behavior of the peak 
height, hx, and the peak width at half maximum, w^- 
Results are shown in Fig. O The peak height/width in- 
creases/decreases with N, both dependences being ap- 
proximately of the power-law type for very large A'^. The 
width decrease is faster than the height increase, so the 
product hxWx decreases. Since the product approximates 
the peak area, its A'^ — > 00 limit should coincide with 
Qx(Xx)- We may therefore conclude that numerical data 



4 



shown in Fig.[2]are compatible with the theory presented 
above, though the convergence to the asymptotic regime 
(presumably of a power-law type) is extremely slow. Let 
us stress that the same type of behavior is expected in all 
sb bosonic models in transition between their 0(2i+2) 
and U(2L+1) dynamical symmetries [l3l |. 

The above-described classification of QPTs in terms of 
nonhermitian degeneracies appears to be very similar to 
the classification of thermodynamic phase transitions in 
terms of zeros of the partition function Z in a complex- 
extended temperature plane [13, HH . This analogy has 
been proposed in Ref. @, but it remained just a sur- 
mise. Here, we can show that the two classifications are 
in fact identical. Recall basic expression for the free en- 
ergy U = — T In Z, specific heat C — —T-g^U, and latent 
heat Q — J C dT, where T is temperature. If neglect- 
ing unimportant prefactors, one sees that these quanti- 
ties are in the same relations as functions Uk, C^, and 
Qk from Eqs. dS])-®. This hnk is vahd if formally as- 
sociating the partition function Z with an appropriate 
power of the partial discriminant Dk , hence also zeros of 
Z with the A:*'^sheet degeneracies. Therefore, the quanti- 
ties introduced here purely on the basis of the Coulomb 



analogy for nonhermitian degeneracies have direct ther- 
modynamic counterparts used in the description of stan- 
dard phase transitions. For instance, the criterion for 
the "force discontinuity" in the first-order QPT, Qk ^0, 
is just the familiar rule of nonzero latent heat. 

In conclusion, we have presented a method to mea- 
sure the distribution of nonhermitian degeneracies close 
to real values of the control parameter for Hamiltoni- 
ans of the form |T]). It is based on electrostatic intu- 
ition, associating degeneracies in the complex plane with 
point charges, and makes it possible to separate degen- 
eracies on Riemann sheets corresponding to different lev- 
els. We formulated criteria for the first-order and contin- 
uous QPTs affecting an arbitrary level and tested them 
in the interacting boson model. Finally, we explained the 
analogy between nonhermitian degeneracies and zeros of 
partition function. These general results can be applied 
to nuclear, molecular, optical, and mesoscopic systems. 
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